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Abstract 

Incidence of cholera outbreak is a serious issue in underdeveloped and developing countries. In Zimbabwe, after the 
massive outbreak in 2008-09, cholera cases and deaths are reported every year from some provinces. Substantial number of 
reported cholera cases in some provinces during and after the epidemic in 2008-09 indicates a plausible presence of 
seasonality in cholera incidence in those regions. We formulate a compartmental mathematical model with periodic slow- 
fast transmission rate to study such recurrent occurrences and fitted the model to cumulative cholera cases and deaths for 
different provinces of Zimbabwe from the beginning of cholera outbreak in 2008-09 to June 2011. Daily and weekly 
reported cholera incidence data were collected from Zimbabwe epidemiological bulletin, Zimbabwe Daily cholera updates 
and Office for the Coordination of Humanitarian Affairs Zimbabwe (OCHA, Zimbabwe). For each province, the basic 
reproduction number {R 0 ) in periodic environment is estimated. To the best of our knowledge, this is probably a pioneering 
attempt to estimate Rq in periodic environment using real-life data set of cholera epidemic for Zimbabwe. Our estimates of 
Ro agree with the previous estimate for some provinces but differ significantly for Bulawayo, Mashonaland West, 
Manicaland, Matabeleland South and Matabeleland North. Seasonal trend in cholera incidence is observed in Harare, 
Mashonaland West, Mashonaland East, Manicaland and Matabeleland South. Our result suggests that, slow transmission is a 
dominating factor for cholera transmission in most of these provinces. Our model projects 6340(5565 — 7264) cholera cases 
and 271(238-309) cholera deaths during the end of the epidemic in 2008-09 to January 1, 2012. We also determine an 
optimal cost-effective control strategy among the four government undertaken interventions namely promoting hand- 
hygiene & clean water distribution, vaccination, treatment and sanitation for each province. 
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Introduction 

Cholera is still a burning problem in underdeveloped and 
developing countries causing morbidity and mortality. In Zim- 
babwe, one of the most severe cholera outbreaks occurred in 
2008-2009, that had been attributed as the worst African 
outbreaks in terms of its high case fatality rate (CFR) and short- 
time extensive spread in some provinces. The outbreak, beginning 
in Chitungwiza, had duration from August 2008 to July 2009, 
ultimately ended with 98,592 reported cases and 4,288 reported 
deaths [1]. These massive outbreaks happened mainly due to 
Zimbabwe's poor health care system, shortage of good-quality 
food and clean drinking water [2] . An economic crisis within this 
period accelerated the deterioration of the country's infrastructure, 
including a breakdown of basic municipal services (such as sewage 
treatment and water supply in many areas) and medical facilities 
[3]. 

The provinces of Zimbabwe experienced a total of 2 1 0 1 cholera 
cases over the period, 17th October, 2009 to 30th June, 2011 
[4,5] . The substantial number of cholera cases in some provinces, 
e.g. Manicaland, Mashonaland West, Masvingo, Midlands, etc., 
both during and after the epidemic in 2008-09, indicate a 



plausible presence of seasonal forcing in cholera incidence in some 
of the provinces. 

Well strategic deployment of cholera intervention/interventions 
in Zimbabwe may reduce future cases and deaths, although the 
projected effect of available cholera interventions is debatable [6]. 
A lot of suggestions have come out for preventing the cholera 
outbreak in those regions. Many regional and international 
organizations suggest providing clean water, hand-hygiene (Soap) 
promotion and construction & promotion of sanitary systems. 
Other groups are arguing for the vaccination program, although 
some experts suggest that the effect of vaccination will be modest 
[7]. Several professionals have also recommended usage of 
rehydration therapy for mild infections (< 10% bodyweight loss) 
and usage of antibiotics (Erythromycin, Doxycycline and Ringer 
Lactate) for severe cases (> 10% bodyweight loss) to reduce 
morbidity [8], cost of productive time loss due to illness, and 
bacterial shedding [9] . With proper treatment of cholera cases, the 
CFR should remain below 1% [10]. However, in terms of cost 
effectiveness, cholera vaccination is by far most costiy interven- 
tion[l 1] with US$1,658 to US$8,274 yields one DALY (Disability- 
adjusted life year) and gaining that same year through promoting 
hand-hygiene need US$3.35, making hand-hygiene the cheapest 
among cholera interventions [12]. Even though, promoting hand- 
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hygiene heavily depends upon the availability of clean water. So 
an optimal balance among different types of interventions may 
significandy reduce the number of cholera cases and deaths at a 
minimal cost. Thus, a well-coordinated effort and an effective 
response to control an outbreak are the most important tasks. 

To control future epidemics, a good understanding of cholera 
transmission dynamics is crucial and mathematical models can be 
utilized as a potential tool [6,13-15]. Some earlier studies on 
cholera are based on the assumption of constant transmission rate 
between human and bacterial population over time [6,13-15] but 
in food or waterborne infections, the role played by temporal 
forcing is more subtle and interesting. There is strong evidence 
that the multi annual dynamics of cholera are interlinked with 
long-term environmental factors [16-18]. 

To capture the presence of possible seasonal pattern within the 
data of reported cholera cases, we include the periodicity factor in 
our model. With these backdrops, we modified the model 
proposed by Hardey el al. [15] to include periodic slow-fast 
transmission and fitted to Zimbabwe's weekly cholera seasonal 
data starting from 2008-2009 epidemics to June 2011. Daily and 
weekly data were collected from Zimbabwe epidemiological 
bulletin [4], Zimbabwe Daily cholera updates [5] and Office for 
the Coordination of Humanitarian Affairs Zimbabwe [19]. The 
basic reproduction number (Rq) carries information about the 
persistence of a disease [20,21]. It is inversely proportional to the 
mean age of (first) infection; greater it is shorter the generation 
time, and the disease transmission will be more explosive [21,22]. 
Aforesaid data was used to estimate Rq in periodic environment, 
for all provinces across the country. To the best of our knowledge, 
this is probably the pioneering attempt to estimate Ro in periodic 
environment using the real data set of cholera epidemic. We 
perform a statistical test suggested by Roger [2 3], using weekly 
cholera incidence data from each province to justify the presence 
of seasonal trend. We also study the existence of any underlying 
pattern of temporal forcing in slow-fast transmission rate with the 
seasonality in cholera incidence, as observed in some provinces. 
We provide forecasts of cumulative cases and deaths from the end 
of epidemic in 2008—09 to 1'" Jan 2012 for different provinces in 
Zimbabwe and study the optimal intervention strategy/strategies 
by minimizing the cost of different cholera interventions. 

Materials and Methods 

Basic model structure 

We modify the existing model [15] assuming temporal 
variations in two types of transmission rates (slow and fast). The 
existing model [15] assumes constant human population size (birth 
and death rates are equal), neglects the cholera-related death rate 
and assumes life time natural immunity to cholera if recovered. 
We have modified these assumptions in our model by incorpo- 
rating variable human population size, cholera-related death rate 
and the effect of natural immunity loss to cholera (as it is now 
proven fact that natural immunity to cholera varies from less than 
one year to two years [24]). Our basic model is a system of five 
differential equations (see Equation 1 & 2) describing how 
individuals can move between different states of susceptibility or 
infection with cholera. 

We categorize the total human populations at time t (denoted 
by N(t)), into susceptible (S(t)), infected (1(f)), and recovered 
(R(t)) classes. A constant recruitment rate (JIh) to human 
population, which is the product of human birth rate (fij,) and 
initial entire human population size (JV(0)), is considered. 



Individuals die naturally at a rate fij. All newly recruited 
individuals are assumed to be susceptible. 

A Recent study showed that freshly shed V. cholerae from human 
intestines are short-lived and hyper-infectious in nature [25] . It out 
competes other V. cholerae grown in vitro, by as much as 700-fold 
for at least the first 5 to 18 hours in the environment [25]. After 
the hyper-infectious stage, V. cholerae organisms lose their 
competitive advantage and become low-infectious. This hyper- 
infectivity is a key factor to understand the explosive nature of 
human-to-human transmission in cholera outbreaks. Based on this 
fact, we classify the bacterial populations in two states, one hyper- 
infectious (Bh) state for fast transmission and this hyper-infectious 
bacteria decay to become low-infectious (B£) state after some 
time, which causes slow environmental transmission. 

Susceptible individuals gain infection by consuming water 
contaminated with the low-infectious bacteria (Bi) and the 



high-infectious bacteria (Bh) at rates Xl - 



t H (t)B H 



K L + B L 



and 



, respectively. The subscripts L and H denote 



Kh + B h 

low infectious and high infectious cholera transmission. Here, 
and Kj] are half saturation constants of low-infectious and high- 
infectious bacterium respectively [15]. Pl(0 = Plo{^ + ^ 



cos(^))} are the 



(1+C0S(^))} and M0 = W1 + ^(1- 

rates of ingesting low-infectious and high-infectious V. cholerae 
bacterium from the contaminated water, which are assumed to be 
time periodic with period 52 weeks. f! L0 , fi m denote the minimum 
transmission rate of low and high infectious V. cholerae respectively 
from the contaminated water, and <5 denotes the amplitude of 
seasonality. Infected individuals either have a natural death (at a 
rate fij) or die due to extreme loss of fluid from their body during 
infection (at a rate fi c ) or recover naturally from cholera infection 
(at a rate y). Recovered individuals are immune to reinfection, but 
this immunity wanes over time and eventually returns to the 
susceptible stage (at a rate CO). During the period of infection; 
infected individuals excrete V. cholerae into water reservoirs around 
them (at a rate £). Since this bacterium is coming directly from 
infected human intestines, it is in hyper-infectious stage and 
decays, ultimately leads to low infectious stage (at a rate x)- The 
natural death rate of low infectious bacterium is &l- 

Based on the above-mentioned assumptions, we construct the 
following system of non-linear differential equations: 



dS 
dt 

dl 
dt 

dR 
dt 



dB L 
dt 



Kh + Bh 
= yI-(a+n d )R 

= U-yB H 

= xBh-SlB l 



K L + B L 



(1) 



where, 
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MO = fco{l + ^(l+cos(^))} 



Pl® = fe 0 {l + |(l+cos(^))} 



(2) 



Model parameters and their interpretations with some parameters' 
base values, taken from previous studies, are given in Table S 1 . A 
flow diagram of the basic cholera model (1) is also given in 
Figure 1. 

Source of data 

Daily and weekly cholera incidence data for each province of 
Zimbabwe have been collected from Zimbabwe epidemiological 
bulletin [4], Zimbabwe Daily cholera updates [5] and Office for 
the Coordination of Humanitarian Affairs Zimbabwe [19] starting 
from 2008-2009 epidemics to June 2011. The data at the 
beginning of the epidemic are quite noisy. To smoothen out the 
initial fluctuations, the data is converted to weekly reported cases 
by aggregating daily case reports for the entire duration of the 
epidemic. Table 1 contains information about the starting points 
and the end points of the data for each province. Total numbers of 
data points vary across the provinces. 

Model calibration 

To calibrate the basic model (1), we have considered the weekly 
reported cholera cases and deaths from each province of 
Zimbabwe starting from 2008-2009 epidemics to June 2011. 
The model is fitted to the cumulative number of cases and deaths 
obtained from the weekly counts in each province. 



The key parameters estimated from the data are the average 
transmission rate of hyper-infectious Bacterium (fino), the average 
transmission rate of low-infectious Bacterium (p L0 ), the amplitude 
of seasonality (<5), the mortality rate of human due to cholera 
infection (fi c ) and the excretion rate of cholera infected individual 
(£). It is not realistic to assume the entire population of a province 
to be susceptible to cholera, as outbreaks generally occur in that 
part of the province where the basic amenities like proper drainage 
system, clean water and food are lacking. So, we first estimate the 
initial number of susceptible, infected and recovered human 
populations from the data by bounding the initial total human 
population size {N(0)) by the total population size of the province. 
The initial concentrations of hyper-infectious (Bh(0)) and low- 
infectious (Bl(0)) bacterium are estimated from the data. Our 
work also involves estimation of initial reported cases (C(0)) and 
deaths (-D(O)) from the data since in some provinces the exact 
reported cases and deaths from the beginning of the epidemic were 
unknown due to reporting delays. 

The cumulative cases and cumulative deaths from the cholera 
model (1) are given by: 



C(t,0) = C(0) + 
D(t~0) = Z>(0) + /i ( , 



S(s)(A L (s) + l H (s))ds 



I(s)ds 



(3) 



where 9eR contains all the unknown variables of the model (1). 
We have n observations (cumulative cases and cumulative deaths) 
from our data at n different weeks as Y(tj) = (C(tj),D(tj)) T , 
where i = \,2,...,n and ti is the i ,h week in our data. 



tta 
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(I 



I 
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Figure 1. Cholera transmission model without any interventions. 

doi:1 0.1 371 /journal.pone.0081 231 .gOOl 
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Table 1. Data Summary. 





Province 


Start date 


End date 


Number of data points 


References 


Harare 


August 18, 2008 


September 12, 2010 


48 


First 12 points from [19], next 30 
points are from [5] and remaining 
6 points are from [4] 


Bulawayo 


November 14, 2008 


April 25, 2009 


17 


First 6 points from [19] and 
remaining 1 1 points are from [5] 


Mashonaland West 


September 21, 2008 


March 27, 201 1 


61 


First 1 1 points from [19], next 29 
points are from [5] and remaining 
21 points are from [4] 


Mashonaland Central 


November 14, 2008 


May 30, 2010 


34 


First 5 points from [19], next 24 
points are from [5] and remaining 
5 points are from [4] 


Mashonaland East 


October 6, 2008 


March 13, 2011 


36 


First 8 points from [19], next 25 
points are from [5] and remaining 
3 points are from [4] 


Midlands 


November 11, 2008 


January 23, 2011 


39 


First 8 points from [19], next 24 
points are from [5] and remaining 
7 points are from [4] 


Masvingo 


November 13, 2008 


June 26, 201 1 


55 


First 5 points from [19], next 33 
points are from [5] ctnd remaining 
17 points are from [4] 


Manicaland 


November 1, 2008 


June 12, 2011 


85 


First 9 points from [19], next 34 
points are from [5] and remaining 
42 points are from [4] 


Matabeleland South 


November 13, 2008 


April 4, 2010 


23 


First 8 points from [19], next 12 
points are from [5] and remaining 
3 points are from [4] 


Matabeleland North 


December 25, 2008 


June 20, 2009 


13 


[5] 



doi:1 0.1 371 /journal.pone.0081 231 .t001 

We assume independent Gaussian prior specifications for 9: 
Oj~N(Yj,r,]), j=\X-4 



(4) 



Let s be the error when fitting cumulative quantities 
(C(t,6),D(t,9)) from the model (1) to the observed data. Then e 
follows independent Gaussian distribution having unknown 
variance o- 2 i.e. 8~N(0,Ia 2 ). For the error variance a Gamma 
distribution is used as a prior for its inverse: 

P (G- 2 )~r('- 



-2\_ rv"° n °-S 2 ) 



2 2 



(5) 



where the prior parameters S 0 and «o m (5) can be interpreted as 
the prior mean for a 1 and the prior accuracy as imaginary 
observations. 

We construct the sum of squares function as: 



estimate the unknown variable 8 for the model (1). Geweke's Z- 
scores [29] were examined to ensure the chain convergence. 

The advantage of using the cumulative over the weekly number of 
new cases in model calibration is that the former smoothes out known 
reporting delays on weekends and national holidays [14,30]. 

Seasonality 

To justify whether the existence of any kind of seasonal forcing 
influence the number of cholera incidence in Zimbabwe provinces, 
a suitable statistical testing procedure is very much needed. The 
weekly data from 18th August, 2008 to 30th June, 2011 is taken 
into account for this purpose. We follow the test procedure 
suggested by Roger [2 3]. 

The entire span of almost 3 years is divided into 52 classes, 
corresponding to the 52 weeks of an year (1 week starting from 
18''' August, 2008) and the total number of cholera cases in week i 
(i= 1,2,. .,52) is denoted by Nj. The probability that any one event 
belongs to i >h class is Pi, where 

qi(\+asi + flci) 



SS(&) = ^[(Citd-CitiM + md-D&M 2 



(6) 



E 4,(1 + as, + 



T,2,..,52. 



(7) 



Posterior distribution of the model unknown variable 6 is 
generated using Delayed Rejection Adaptive Metropolis algorithm 
(DRAM) [26,27] with an initial burn of 100000 iterations. MCMC 
toolbox in MATLAB written by Marko Laine [28] was used to 



where, qi denotes the frequency for class i under the null hypothesis, 
Sj= sin (2ni/52), c, = cos (2ni/ 52), a. and fi are the parameters of 
the model (7). Hq : ot = /> = 0 indicates the absence of seasonality and 
Hi : txj^O or [1^0 indicate the seasonality in cholera incidence. 
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The test statistics for testing Hq is of the form Let, Y(t,s), t>s be the evolution operator of the linear w- 

periodic system 

A 2 u+B 2 v-2ABw 

n(uv-w 2 ) ■ («) * = _ v{t)y (io) 



Where, 



That is, for each seM., the 3x3 matrix Y(t,s) satisfies 



A = ^2 DiSi =^2(Ni- nqi)s h 

i i 

B=J2 Did =J2(Ni- nqi)a, 

i i 

w = Yl qiS,Ci ~ (y! qiS ^j (y! qiC ^j ' 

and n = y\iV,-. 

The test statistic R is asymptotically distributed as chi-square 
with 2 degrees of freedom. 

Estimating reproductive numbers in periodic 
environments 

The Model (1) has a unique disease free equilibrium given by: 



E 0 = (5^,0,0,0,0). 



(9) 



Following [31], we calculate the matrix of new infection from 
our system (1) as: 



F(t)- 



0 



0 0 0 



0 0 



0 



and the transmission matrix as: 



V(t)- 



(y+H c + n d ) 0 0 

-e x o 

0 -x Sl 



dt 



Y(t,s) 



-V(t)Y(t,s), 



for all t~>s and Y(s,s) = I, where I is the 3x3 identity matrix. 

Let Cnj be the ordered Banach space of all K7-periodic functions 
from R to IR 3 which is equipped with maximum norm and 
the positive cone = {<j)eC w : (j)(t)>0, for all t in IR}. Consider 
the following linear operator L : C ro -> C w by 



(I4)(t) 



J+ 00 Y(t, t - a)F(t - a)(j>(t -a)da. (11) 



Following Wang and Zhao (2008) [31], we call L the next 
infection operator, and define the basic reproduction number (Rq) 



p{L) 



(12) 



where p{L) is the spectral radius of the operator L defined in 
Equation (11). 

Motivated by the concept of the partial reproduction numbers 
defined by Mukandavire et.al. [14], we similarly define two partial 
reproduction numbers Ri and R/ t in periodic environment. The 
subscripts / and h correspond to low infectious and high infectious 
transmission, respectively. 

Basic reproduction number (Rq) is the sum of two partial 
reproductive numbers- the one is arising from the contact between 
the human and low-infectious bacteria, which we denote as R/ and 
the other one arising from the contact between the human and 
hyper- infectious bacteria, denoted as R/ t . Using Lemma 1 given 
in Appendix SI and the estimated parameter values (Table S3 
and Table S4), we numerically estimate Rq, Ri and Rh for each 
province. Details of the derivation procedure oi Rq, Ri and R/, are 
given in Appendix SI. For uncertainty, we draw 95% confidence 
interval around the estimated values. The following procedure is 
applied to derive the 95% confidence intervals for Rq, Ri and Rh, 
respectively. 

We draw a sample of size n (n=\ 000) from the posterior 
distribution of 9 (set of model variables, which are estimated) using 
simple random sampling without replacement (SRSWOR) 
scheme. The posterior distribution of 6 is depicted in Figure SI 
and Figure S2. For each of the sample values of 9, we estimate 
numerically the value of Rq (using Lemma 1, Appendix SI). 
Thus, a vector of size n, is generated for Rq. Curtailing the lower 
2.5% and the upper 97.5% observations from the ordered vector 
of Rq, we obtain the 95% confidence interval for Rq. Applying the 
similar procedure we draw 95% confidence intervals for two 
partial reproductive numbers, Ri and Rh, respectively. 

Projection of future cases and deaths 

We project the number of cholera cases and deaths from the 
end of epidemic in 2008-09 to January 1, 2012. For uncertainty, 
we derive 95% credible intervals around the estimates of future 



PLOS ONE | www.plosone.org 



5 



December 2013 | Volume 8 | Issue 12 | e81231 



A Study on Cholera Seasonal Data of Zimbabwe 



projected cases and deaths. To predict the number of cases and 
deaths for a particular province, we simulate the cholera model (1), 
using the known & estimated parameters (Table SI and Table 
S3) and demographic parameters (Table S4), up to the end of 
epidemic in 2008-09, in that region. We obtain different 
demographic variables of human & pathogen (S, I, R, Bfj and 
B£), new cases and deaths corresponding to the end of epidemic in 
2008-09. Using this information from the previous simulation as 
initial conditions and parameter values from Table SI and 
Table S3, we simulate the model (1) to obtain predicted cases and 
deaths from the end of the epidemic in 2008-09 to January 1, 
2012. 

We used the following procedure to derive the 95% confidence 
intervals for projected cases and deaths. For each of the sample 
value of 9 (see, section-Estimating reproductive numbers in 
periodic environments), we predict the number of cases and 
deaths using the above procedure. Thus, two vectors, each of size 
n, axe generated for predicted cases and deaths, respectively. 
Curtailing the lower 2.5% and the upper 97.5% observations from 
the ordered vector of predicted cases and deaths we obtain the 
95% confidence intervals for projected cases and deaths, 
respectively. 

Model with different cholera interventions 

Effect of four different types of cholera interventions namely 
hand-hygiene promotion & clean water supply, treatment using 
oral rehydration therapy & antibiotics, vaccination and sanitation 
are studied. We assume that hand-hygiene (soap) & clean water 
will reduce bacterial ingestion by a fraction (1 — 9{t)), where 9(f) is 
the relative rate of reduction in bacterial ingestion per week using 
hand-hygiene & clean water supply. Vaccinated population is 
increased by a proportion p(t)o of the susceptible individuals, who 
are successfully vaccinated, where p(t) is the per week vaccination 
rate and a is the vaccine efficiency. Vaccinated population is 
decreased due to the waning of vaccine based immunity (at a rate 
fi) to become susceptible again and die (natural deaths) at a rate fi d . 
We assume that a proportion a(t) of the infected individuals 
receive treatment by oral rehydration salt (for < 10% body weight 
loss) and by antibiotic (for > 10% body weight loss) per week. 
Natural recovery rate of the treated person increases by relative 
rate of recovery X. Since, excretion can be affected by the use of 
antibiotics [32,33], the relative rate of shedding is reduced by a 
fraction \\i among the proportion of the infected individuals who 
receive antibiotic at a rate a(t) per week. Now proper sanitary and 
drainage system will prevent the human waste to contaminate the 
nearby water reservoirs. Invariably, sanitation will reduce excre- 
tion rate of human that contaminate the nearby reservoirs by a 
fraction (1 — s(t)), where s(t) is the rate of reduction in human 
shedding per week by construction and promotion of sanitation. In 
Zimbabwe, 80% of the total populations have access to improved 
water source and 40% of the total populations have access to 
proper sanitation facilities [1]. Therefore, we assume maximum 
percentage reduction in bacterial ingestion rate through hand 
hygiene & clean water supply and reduction in human shedding 
possible by promoting sanitation in a week to be 80% and 40%, 
respectively. We also assume that maximum 70% of the infected 
individuals receive proper treatment in a week and maximum 
vaccination coverage possible in a week is about 35% of the total 
susceptible population [34]. Effects of hand-hygiene & clean water, 
vaccination, treatment, sanitation and their different combinations 
are projected from the end of epidemic in 2008-09 to January 1, 
2012. 

System of non-linear differential equations representing the 
effect of different interventions on our basic model (1) is given as 



follows: 

dS 
dt 



= n H +wR+ev-(i-0(t))[p H (t) 



B H S 
Kh + B h 



d£ 
dt 



dV 
dt 

dR 
dt 

dB H 
dt 

dB L 
dt 



+ P L (t)-^^\-p(t)aS-!i d S 

- G«rf + n c + 0 - <t))y + <t)yW 

= ptfaS-eV-fijV 

= ((l-x(t))y + x(t)yA)I-(oj + p d )R 

= (l - s (t))(Mt)+(i-<t)Mi-xB H 

= xB H -5 L B L 



(13) 



Intervention parameters and their interpretations with some 
parameters' base values taken from earlier studies are given in 
Table S2. A flow diagram of the intervention model is depicted in 
Figure 2. 

An optimal intervention strategy 

To determine the optimal intervention strategy/strategies 
(which reduce the number of cases and deaths projected from 
the end of the epidemic in 2008-09 to January 1, 2012, at a 
minimal cost), we define the following cost function: 

J(9,p,a,s) = \ T {An c I(t) + {Bp(t)S(t) + Cp 2 (t)<j} 
Jo 

+ {Doc(t)I(t) + Ea 2 (t)} ( 14 ) 

+ {F9(t) + G8 2 (t)\ + {Hs(t) + Ks 2 (t)}]dt 

where interventions are applied for T weeks. First term in the 
right-hand side of (14) represents the cost of cholera-related deaths 
and remaining terms are costs associated with the implementation 
of different interventions. Nonlinear terms in the objective 
function / represent the costs of interventions in emergency 
situations. A, B, C, D, E, F, G, H and K are fixed cost 
coefficients, given in the Table 2. 

Fixed costs are transformed to subsequent intervention year 

Y 

costs by multiplying with a constant c, where c = — . 

* cost year 

Y intervention and Y cost year are the annual consumer price indices 
(US) when interventions are applied (i.e. from the end of epidemic 
in 2008-09 to January 1, 2012) and the annual consumer price 
index (US) of the years (see Table 2) respectively. Y interw „ tion and 
Y cost year are in the same base 1982=100, data of annual 
consumer price index (US) were collected from U.S. Bureau of 
Labor Statistics [35]. 
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Figure 2. Cholera transmission model with different interventions. 
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Our goal is to minimize the objective function J with respect to 
different control parameters 8(t), p(t), oc(t) and s(t) to determine 
an optimal intervention combination. This is a dynamic control 
problem and is solved directly by using the Pontryagin's 
Maximum Principle [36] and the method of steepest decent 
[37]. Minimization procedure of the objective function / is briefly 
described in Appendix SI. 

The average coverage percentages (per week) of hand-hygiene 
(soap) & clean water distributions, treatment and sanitation are 
estimated using the following formulas: 



dt 



x 100, oc= 



dt 



x 100 and 



dt 



x 100, 



where, 6, a. and s, denote the average coverage percentages (per 
week) of hand-hygiene (soap) & clean water distributions, 
treatment and sanitation, respectively. 9 c (t), a c (?) and s c (t) are 
the optimal rates of hand-hygiene (soap) & clean water distribu- 



Table 2. Fixed cost-coefficients. 




Notations 


Interpretations 


Year 


Value (US$) 


Reference 


A 


Cost of productive time lost per premature death (Calculated 
with life expectancy 73 years) 


2006 


27264.6 


[56] 


B 


Cost of oral cholera vaccine (OCV) per fully immunized person 


2004 


7 


[57] 


C 


Cost of OCV per fully immunized person in high emergencies 


2005 


18 


[57] 


D 


Cost of medicines and health centre consultation per mild/moderate case 


2006 


23.7 


[56] 


E 


Cost of medicines and hospital admission per severe cholera cases 


2006 


70.2 


[56] 


F 


Cost of per percent reduction in bacterial ingestion rate by promoting 
hand-hygiene and water supply 


1993 


320.35 


[12] 


G 


Cost of per percent reduction in bacterial ingestion rate by promoting 
hand-hygiene and water supply in high emergencies 


1993 


448.49 


assumed 40% increase in 
normal cost 


H 


Cost of per percent reduction in human shedding by promoting sanitation 
(construction and promotion of latrine and drainage system) 


1993 


270 


[12] 


K 


Cost of per percent reduction in human shedding by promoting sanitation 
(construction and promotion of latrine and drainage system) in high emergencies 


1993 


378 


assumed 40% increase in 
normal cost 


doi:1 0.1 371 /journal.pone.0081 231 .t002 



PLOS ONE | www.plosone.org 



December 2013 | Volume 8 | Issue 12 | e81231 



A Study on Cholera Seasonal Data of Zimbabwe 




1000 



x 10 



60 120 



10000 



60 



140 



10000 



60 



140 



xlO 



60 



140 



; 10000 



40 

Time (weeks) 



80 




2000 



Time (weeks) 



Figure 3. Province-wise cumulative cholera cases in Zimbabwe. 

The observed data points (available at some discrete time points over a 
time period, which varies across the study regions) are shown by blue 
circles while the solid lines depict the model solutions. The cumulative 
cholera cases from the model are plotted for each day of the time 
period (from the start to end week for the observed cholera data) using 
parameter values and initial conditions from Table S3 and Table 54. 
The above plots of cholera cases from the different provinces of 
Zimbabwe are as follows: (i) Harare; (ii) Bulawayo; (iii) Mashonland West; 
(iv) Mashonland Central; (v) Mashonland East; (vi) Midlands; (vii) 
Masvingo; (viii) Manicaland; (ix) Matabalend South; and (x) Matabalend 
North. 

doi:1 0.1 371 /journal.pone.0081 231 .g003 

tions, treatment and sanitation, respectively, for which the cost 
function J (see equation (14)) is minimum. T denotes the total 
number of weeks during which an intervention is applied. 

Total coverage percentage of vaccination is estimated using the 
following formula: 

t-T 



1/(0) + 



Total vaccination coverage 



p c (t)aS{t)dt 
£ x 100, 



S(t)dt 



where, V(0) is the initial number of vaccinated individuals and 
S(t) is the number of susceptible individuals at week /. a is the 
vaccine efficiency (Table S2). p c (i) is optimal vaccination rate 
which minimizes the cost function / (see equation (14). 

Cost per averted case for an intervention is calculated using the 
following formula: 



Cost per averted case 

Total cost of an intervention 
Total number of cases averted by that intervention 




40 

Time (weeks) 



Time (weeks) 



Figure 4. Province-wise cumulative cholera-related deaths. The 

data points are shown by empty blue circles while the model fits by the 
solid lines. The plots are given in the same order as of Figure 3. The 
cumulative deaths from the model are plotted using parameter values 
and initial conditions from Table S3 and Table S4. 
doi:1 0.1 371/joumal.pone.0081 231 .g004 

A 95% confidence interval for each of the following quantities 
are obtained following the same technique as explained in sections 
Projection of future cases and deaths and Estimating 
reproductive numbers in periodic environments: (1) cases 
that occurred in spite of applying an intervention, (2) total cost of 
an intervention and (3) cost per averted case. 

Results 

Cholera model fitting for the cumulative reported cholera cases 
and deaths are depicted in Figure 3 and Figure 4, respectively. 
A comparison between weekly reported cholera cases & deaths 
from each province with the model solution are shown in 
Figure 5(A, B, and C) and Figure 6(A, B, and C), 
respectively. The estimated model parameters, including human 
and pathogen demographic parameters, for each province are 
given in the format [estimate (95% CI)], (Table S3 and 
Table S4). Plots for the posterior distributions of the estimated 
unknown variables of the cholera model (1) are given in Figure 
SI and Figure S2. 

The contributions of low-infectious {Ri) versus high-infectious 
(R/i) transmission to Rq vary widely. Our estimated values of Rq, 
R/ and R/, are in good agreement with the previous estimates 
given by Mukandavire et. al. [14], for the provinces Harare, 
Mashonaland East, Mashonaland Central, Midlands and Mas- 
vingo but significantly differ in Bulawayo, Mashonaland West, 
Manicaland, Matabeleland South and Matabeleland North. In 
Mashonaland West and Manicaland, contribution of low-infec- 
tious (Ri) is higher than hyper-infectious (R/,) transmission to Rq 
(see Table 3). Opposite trend is observed in the estimates of Ri and 
R/i given by Mukandavire et. al. [14] in these two provinces. 
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Figure 5. Cholera model fitting for the weekly new cholera cases. The solid line represents the model solution, and blue circles mark the 
reported cholera cases in the provinces using parameter values and initial conditions from Table S3 and Table S4. 
doi:1 0.1 371 /journal.pone.0081 231 .g005 
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Figure 6. Cholera model fitting for the weekly new cholera deaths. The solid line represents the model solution, and blue circles mark the 
reported cholera deaths in the provinces using parameter values and initial conditions from Table S3 and Table S4. 
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Table 3. Estimates of R/, R/, and R 0 . 





Zimbabwe province 


Ri 


95% CI 


% Ro 


Rh 


95% CI 


% Ro 


Ro 


95% CI 


Harare 


1.2294 


1.1885-1.2764 


99.06 


0.011 


0.0004-0.04 


0.94 


1 .2406 


1 .206-1 .285 


Bulawayo 


0.098 


0.0212-0.2083 


95.89 


0.0042 


0-0.025 


4.11 


0.1022 


0.0215-0.2382 


Mashonaland West 


1.4105 


1.338-1.4651 


97.66 


0.0338 


0.0017-0.1387 


2.34 


1.4443 


1.4126-1.4782 


Mashonaland Central 


0.005 


0.0044-0.0057 


0.28 


1.7701 


1.7465-1.8015 


99.72 


1.775 


1.7518-1.8063 


Mashonaland East 


0.0003 


0-0.0015 


0.0163 


1 .8453 


1.8387-1.8499 


99.98 


1 .8456 


1.8397-1.8501 


Midlands 


0.0071 


0.0059-0.0092 


0.39 


1 .7974 


1.7844-1.8101 


99.61 


1 .8045 


1.7929-1.8161 


Manicaland 


1.1146 


1.031-1.1575 


97.91 


0.0238 


0.0012-0.1256 


2.09 


1.1384 


1.0982-1.1765 


Masvingo 


0.0012 


1 .52E-04 - 0.003 


0.065 


1.8235 


1.8139-1.8350 


99.94 


1 .8246 


1.8158-1.8355 


Matabeleland North 


0.0353 


0.0164-0.0875 


65.25 


0.0188 


7.29E-04 - 0.0891 


34.75 


0.0541 


0.0208-0.138 


Matabeleland South 


0.4925 


0.4697-0.5158 


62.23 


0.2989 


0.2499-0.3756 


37.77 


0.7914 


0.74-0.8684 
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Estimates of Ro in Bulawayo (0.1022), Matabeleland South 
(0.7914) and Matabeleland North (0.0541) are all found to be 
below the unity, which are drastically different from the estimates 
given by Mukandavire et. al. [14]. 

To justify the existence of any seasonal trends in cholera 
incidence data of different provinces, the values of the test statistic 
R, defined in (8), are calculated using weekly data from 18th 
August, 2008 to 30th June, 2011. The values of R and 
corresponding p-values for each province are given in Table 4. 
Significant seasonal trends in cholera incidence data are observed 
in Harare, Mashonaland West, Mashonaland East, Manicaland 
and Matabeleland South. Among these five provinces, Harare and 
Manicaland exhibit highly significant seasonal trend, as confirmed 
by the corresponding p-value of the test (< 0.01). 

Our basic model (1) without any interventions, projects 
6340(5565-7264) cases and 271(238-309) deaths due to 
cholera in Zimbabwe between the end of epidemic in 2008-09 
to January 1, 2012. Table 5 and Table 6 contain the predicted 
total cholera cases and deaths for the provinces during that time 
interval. Among the ten provinces, in Mashonaland West and 
Mashonaland Central the most numbers of cholera cases and 
deaths are predicted. In Mashonaland West, total 

Table 4. Table showing results for seasonality testing. 



Province 


Test statistic R 


p-value 


Harare 


9.7992 


0.0074" 


Bulawayo 


0.3625 


0.8342 


Mashonaland West 


6.9417 


0.0311* 


Mashonaland Central 


0.1040 


0.9493 


Mashonaland East 


6.1335 


0.0466* 


Midlands 


0.7608 


0.6836 


Masvingo 


1.1860 


0.5527 


Manicaland 


13.0533 


0.0015" 


Matabeleland South 


7.1850 


0.0275* 


Matabeleland North 


0.1395 


0.9326 



Bold provinces are where the seasonality test result is found to be positive. * : 
Denote the provinces where seasonality presents at the significance level 0.05 
and * *: denote the provinces where seasonality present at the significance 
level 0.01. 

doi:1 0.1 371/journal.pone.0081 231 .t004 



3 1 1 8(2788 - 3463) cases, and 1 35( 120 - 1 5 1) deaths are predicted 
during the mentioned period. This is about 49% of the total 
predicted cases and about 50% of the total predicted deaths in all 
provinces of Zimbabwe. In Mashonaland Central, total 
987(858-1161) cases, and 33(29-39) deaths are predicted 
during the mentioned period. This is about 16% of the total 
predicted cases and about 12% of the total predicted deaths in all 
provinces of Zimbabwe. 

To justify the predictive performance of our basic model (1), we 
compare the predicted cumulative cases and deaths from the end 
of the epidemic in 2008-2009 to January 1, 2012 with the 
reported cases and deaths' figures. The reported cumulative cases 
and deaths during the aforesaid time period are 2225 and 72, 
respectively [38], which is about 35%(31% — 40%) of the model 
projected cases and about 27%(23% — 30%) of the model 
projected deaths. Significant difference in the actual and predicted 
case and death figures may be attributed to the higher percentage 
of underreporting of cholera cases and deaths [39], that was not 
considered while making these predictions. According to WHO, 
the officially reported cholera cases represent only 5 — 10% of the 
actual number of cases those are occurring annually worldwide 
[39]. 

We found that, in the African region the countries report 
cholera cases more consistently than the other countries under 
WHO [39]. Also Zimbabwe's Integrated Diseases Surveillance & 
Response Technical guidelines list Cholera among the diseases 
that must be reported on a daily basis during epidemics to prevent 
avoidable illness and death [40]. Thus we may expect that the 
percentage of reported cases is higher in Zimbabwe than the 
worldwide statistics of under-reporting, although, we do not have 
the specific figures/numbers from literature depicting the actual 
percentages of under-reporting in Zimbabwe during the end of 
epidemic in 2008-09 to January, 1, 2012. The actual reported 
cases during that period are about 3 1 %— 40% of our model 
predicted cases. Hence this percentage may be considered as an 
estimate of reporting of cholera cases in Zimbabwe, which is much 
greater than the worldwide statistics (5%-10%). 

It is already mentioned that Mashonaland West and Mashona- 
land Central are high-risk provinces in terms of cholera incidence 
between the end of epidemic in 2008-09 and January 1, 2012. 
Therefore, we discuss the results of different interventions and 
their layered combinations for these two provinces only. 

We have arrived at the following conclusions. 

In Mashonaland West, on average 20.18%(19.35%-22.53%) 
(Table 7) the relative reduction of bacterial ingestion (per week) is 
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Table 6. Number of deaths from cholera projected between the end of 2008-09 epidemic to January 1, 2012, by province under 
base case and under each intervention scenario at an optimal rate. 







Harare 


Mash 
Bulawayo west 


Mash 
cen 


Mash 
east 


Midlands 


Masvin 


Manica 


Mata 
south 


Mata 
north 


Total 


Base deaths 


26 (22-31) 


0 (0-0) 


135 (120- 
151) 


33 (29-39) 


10 (9-12) 


8 (6-10) 


9 (8-11) 


38 (34-42) 


6 (5-6) 


6 (5-7) 


271 (238- 
309) 


PH & CWD 


5 (4-6) 




19 (18-21) 


3 (2-3) 


0 (0-0) 


0 (0-0) 


0 (0-0) 


7 (6-8) 


2 (2-2) 


3 (2-3) 


39 (34-43) 


VA 


24 (8-31) 




107 (32-138)19 (5-34) 


9 (0-11) 


7 (0-10) 


8 (0-11) 


33 (12-42) 


5 (3-6) 




212 (60- 


TR 


8 (7-9) 




31 (29-33) 


3 (2-3) 


0 (0-0) 


0 (0-0) 


0 (0-0) 


12 (11-13) 


2 (2-3) 




56 (51-61) 


SN 


16 (14-19) 


* 


68 (62-73) 


5 (4-6) 


0 (0-0) 


0 (0-1) 


0 (0-0) 


25 (22-27) 


5 (4-5) 




119 (106- 
131) 


PH & CWD + TR + VA 


3 (3-4) 




13 (11-14) 


2 (2-3) 


0 (0-0) 


0 (0-0) 


0 (0-0) 


5 (4-5) 


1 (1-D 


2 (2-3) 


26 (23-30) 


PH & CWD + TR + SN 


3 (3-4) 




13 (12-14) 


2 (2-3) 


0 (0-0) 


0 (0-0) 


0 (0-0) 


5 (4-5) 


1 (1-D 


2 (2-3) 


26 (24-30) 


PH & CWD + VA + SN 


5 (3-5) 




17 (14-19) 


3 (2-3) 


0 (0-0) 


0 (0-0) 


0 (0-0) 


7 (5-7) 


2 (1-2) 


3 (2-3) 


37 (27-39) 


VA + TR + SN 


7 (4-8) 


* 


28 (18-30) 


2 (2-3) 


0 (0-0) 


0 (0-0) 


0 (0-0) 


10 (6-12) 


2 (1-2) 


4 (3-4) 


53 (34-59) 


PH & CWD + TR + 
VA + SN 


3 (3-4) 




13 (11-14) 


2 (2-3) 


0 (0-0) 


0 (0-0) 


0 (0-0) 


5 (4-5) 


1 (1-D 


2 (2-3) 


26 (23-30) 



Data are given in the format [mean (95% CI)]. 

Notations in the first column are exactly same as Table 5. * indicate the intervention/intervention combination which do not have any effect on death reduction in a province. 
doi:1 0.1 371/journal.pone.0081 231 .t006 



observed. By promoting hand-hygiene and clean water supply will 
avert 2782(2478-3082) cases (Table 5) and 116(102-130) 
deaths (Table 6). Cost of carrying out this intervention over the 
period (end of epidemic in 2008-09 to January 1, 2012) is about 
281500(261560-324450) (USD) (Table 8). Cost per averted case 
in Mashonaland West using hand-hygiene and clean water supply 
is 101.42(94.04-122.93) (USD) (Table 9), making this interven- 
tion as the cheapest among other single interventions. Again, in 
spite of total vaccination coverage 0.96(0.0002— 13.49) (Table 7), 
the projected cases and deaths occurred in Mashonaland West will 
be 2446(653-3175) and 107(32-138), respectively (Tables 5 
and 6). Cost per averted case in Mashonaland West using 
vaccination is 197220(2478-685370) (USD) (Table 9), making 
this intervention the most costly among other single interventions. 

Among layered interventions in Mashonaland West, hand- 
hygiene & clean water distribution with treatment and sanitation is 
the most cost-effective and will avert 2815(2507 — 3143) cases 
(Table 5) and 122(108-137) deaths (Table 6). Cost per averted 
case using this intervention is 56.14(51.98-59.93) (USD) 
(Table 9). Hand-hygiene & clean water distribution with 
vaccination and sanitation is the most costly intervention in 
Mashonaland West with cost per averted case is 
2250(88.55-8705) (USD) (Table 9). 

For Mashonaland Central, among single interventions, hand- 
hygiene & clean water distribution will avert most numbers of 
cases and deaths. Projected cases and deaths using hand-hygiene & 
clean water distributions are 50(39 — 61) and 3(2 — 3), respectively 
(Tables 5 and 6). Cost per averted case using this intervention is 
40.16(33.48 -45.89)(USD) (Table 9). Treatment is found to be 
the cheapest among other single interventions in Mashonaland 
Central with cost per averted case being 
30.93(24.09- 36. 75)(USD) (Table 9). Treatment of average 
13.91% (13.39%-15.06%) (Table 7) cholera infected individuals 
(per week) will avert 901(791 — 1054) cases (Table 5) and 
30(27 — 36) deaths (Table 6), respectively. 



Among layered interventions in Mashonaland Central, hand- 
hygiene & clean water distribution with vaccination and sanitation 
will avert most numbers of cases and deaths. Projected cases and 
deaths using this layered intervention are 47(37 — 59) and 3(2 — 3), 
respectively (Tables 5 and 6). Cost per averted case using this 
intervention is 1362.89(32.42-20415.19) (USD) (Table 9). In 
terms of cost effectiveness hand-hygiene & clean water distribution 
with treatment and sanitation is found to be the cheapest layered 
intervention with cost per averted case being 23.58(19.13 — 27.40) 
(USD) (Table 9). This intervention combination will avert 
937(819-1100) cases and 31(27-36) deaths, respectively 
(Tables 5 and 6). 

Discussion 

Our analysis suggests that, the routes of cholera transmission 
vary from province to province, that agrees with the findings of 
Mukandavire et. al. [14]. This heterogeneity in transmission 
dynamics may be due to the diverse geographic and climatic 
conditions across the country. A similar pattern in transmission 
dynamics is observed in Harare, Mashonaland West, Manicaland 
and Matabeleland South, where seasonality in cholera incidence 
was observed. In these provinces, slow transmission route is a 
dominating factor (fi L0 > Puq) for cholera transmission (Table 
S3). Earlier studies on cholera suggest that slow transmission route 
is more correlated to the climatic and environmental factors 
[17,18,41-47] and is the main cause for seasonal dynamics of 
cholera [41,48-51]. Unfortunately, due to lack of climatic data of 
Zimbabwe, we are unable to draw any quantitative inference, for 
example, whether inter-annual climatic variation in different 
provinces affects the transmission dynamics of cholera or not. 

Estimate of human shedding rate (£) (Table S3) in Bulawayo 
province differs from other nine provinces (one order of magnitude 
higher than for the other provinces). A possible reason for such 
difference may be due to the fact that the city receives portable 
water supply from five surface dams [52] and constantly suffers 
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Table 7. Average optimal rate at which different intervention should be given between the end of 2008-09 epidemic to January 1, 
2012, for each province. 





Harare 


Mash west 


Mash cen 


Mash east 


Midlands 


PH & CWD 


12.91(12.16-14.01) 


20.18(19.35-22.53) 


7.65(6.73-9.63) 


3.85(2.37-13.66) 


4.26(3.58-7.24) 


VA 


1.17(2.82E-04-14.95) 


0.96(2.96E-04-1 3.49) 


1.61(4.46E-04-17.36) 


0.96(3.8E-04-13.22) 


1.25(2.52E-04-11.74) 


TR 


23.86(22.88-25.33) 


30.05(29.49-30.79) 


13.91(13.39-15.06) 


4.89(3.44-11.76) 


9.31(8.81-10.96) 


SN 


10.33(9.11-11.63) 


16.66(16.19-17.35) 


5.73(5.23-7.13) 


2.72(2.57-3.05) 


3.82(3.38-5.70) 


PH & CWD + TR + VA 


8.16(4.12-9); 20.32(8.53- 
22.38); 0.63 
(2.78E-04-8.11) 


11.51(5.95-12.34); 
25.11(10.87-27.24); 
0.91 (3E-04-1 3.53) 


2.41(1.07-3.17); 

13.14(4.79-14.48); 

0.95(4.39E-04-15.01) 


0.30(0.30-0.31); 
3.72(1.47-6.45); 
1.57(3.8E-04-10.6) 


0.31(0.24-0.47); 

8.85(2.59-13.52); 

0.97(2.52E-04-13.21) 


PH & CWD + TR + SN 


8.40(8.04-8.85); 

21.21(19.58-22.01); 

0.16(0.15-0.17) 


12.05(11.78-12.25); 
26.65(26.09-27.39); 
0.82(0.74-0.91) 


2.49(1.99-3.13); 
14(13.43-15.99); 
0.15(0.15-0.15) 


0.30(0.30-0.30); 
4.59(3.26-6.98); 
0.0055(0-0.056) 


0.32(0.24-0.55); 

10.08(8.58-21.71); 

0.023(0-0.51) 


PH & CWD + VA + SN 


11.34(5.09-13.48); 

1.54(2.82E-04-13.86); 

1.53(0.43-1.94) 


15.09(7.09-19.40); 
2.78(3E-04-13.99); 
3.07(0.88-4.80) 


7.38(2.18-12.47); 

0.75(4.39E-04-14.46); 

2.09(0.44-6.77) 


3.83(0.90-19.34); 

0.73(3.8E-04-11,43); 

2.5(0.26-14.59) 


4(0.86-8.55); 0.81(2.52E- 
04-11.3); 2.31(0.09-6.32) 


VA + TR + SN 


1.96(2.82E-04-15.35); 

20.24(8.53-23.46); 

2.85(0.81-3.55) 


0.196(2.95E-04-4.099); 

27.85(11.2-29.15); 

5.33(1.44-5.63) 


1.23(4.39E-04-14.38); 

12.66(4.76-14.53); 

0.94(0.44-1.26) 


0.45(3.8E-04-7.31); 

4.31(1.66-7.13); 

0.16(0.15-0.20) 


1.46(2.52E-04-13.40); 

8.4(2.6-13.07); 

0.16(0.13-0.19) 


PH & CWD + TR+ VA + SN 


8.25(4.02-8.97); 
20.85(8.45-22.45); 
0.45(2.39E-04-10.93); 
0.16(0.15-0.18) 


11.32(5.93-12.23); 
24.74(10.85-27.27); 
0.84(3E-04-12.03); 
0.75(0.21-0.93) 


241(1.09-3.21); 
13.23(4.85-15.61); 
0.996(4.39E-04-14.36); 
0.15(0.15-0.15) 


0.30(0.30-0.30); 
4.38(1.83-7.6); 
0.51 (3.8E-04-1 2.08); 
0.005(0-0.05) 


0.31(0.24-0.34); 
10.35(8.79-19.02); 2.7E- 
04(2.5E-04-3.04E-04); 
0.002(0-0.039) 




Masvin 


Manica 


Mata south 


Mata north 





PH & CWD 


4.21(2.37-17.14) 


13.33(12.95-13.79) 


7.16(7.02-7.37) 


6.89(6.59-7.31) 


VA 


1.65(4.28E-04-16.61) 


2.7(2.3E-04-18.8) 


1.19(0.001-12.6) 




TR 


4.97(3.54-7.46) 


24(23.3-24.7) 


13(12.72-13.36) 


* 


SN 


2.64(2.46-3.52) 


10.02(9.25-11) 


3.67(3.5-3.87) 


* 


PH & CWD + TR + VA 


0.28(0.25-0.30); 
4.65(1.47-8.93); 
0.91 (4.28E-04-1 1.4) 


8.85(4.79-9.6); 
20.05(9.18-21.94); 
0.49(2.3 E-04-5.66) 


4.33(2.39-5.05); 

11.08(6.17-12.74); 

2.02(0.001-15.01) 


5.45(4.63-5.90); 

11.52(10.02-12.56); 

0.19(7.0E-04-2.37) 


PH & CWD + TR + SN 


0.28(0.25-0.30); 
4.94(3.46-8.14); 5.49E- 
04(0-0.0074) 


9.35(9.11-9.56); 

21.42(20.93-22.05); 

0.22(0.19-0.25) 


4.88(4.68-5.06); 

12.48(12.16-12.8); 

0.024(0.008-0.0415) 


5.49(5.19-5.96); 

11.6(11.07-12.6); 

5.74(0-7.21) 


PH & CWD + VA + SN 


4.24(0.498-17.86); 

1.22(4.28E-04-11.6); 

2.70(0.48-11.85) 


12.18(5.75-13.34); 

0.83(2.3E-04-11.75); 

2.0(0.53-2.39) 


6.34(3.24-7.3); 

1.59(0.001-13.42); 

0.20(0.11-0.27) 


6.87(6.58-7.38); 9.6E- 

04(6.5E-04-1.1E-03); 

4.14(0-5.18) 


VA + TR + SN 


1.01(4.28E-04-1 1.5); 

4.51(1.46-6.42); 

0.148(0.148-0.149) 


1.51(2.3E-04-15.14); 

21.09(9.32-23.47); 

3.32(0.99-3.79) 


I. 37(0.001-13.72); 

II. 92(6.22-13.14); 
1.02(0.25-1.3) 


0.036(7E-04-0.44); 

11.57(10.06-12.57); 

0.0018(0-0.033) 


PH & CWD + TR + VA + SN 


0.28(0.25-0.30); 
5.09(1.51-10.54); 
0.51 (4.28E-04-1 2.07); 
6.37E-04(0-0.01) 


8.79(4.70-9.54); 
19.91(9.13-21.9); 
0.87(0.0002-12.57); 
0.21(0.15-0.25) 


4.88(4.67-5.06); 
12.47(12.17-12.79); 
0.001(0.001-0.001); 
0.023(0.008-0.0409) 


5.4(2.9-6.0); 1 1 .4(6.7- 
12.6); 0.53(6E-04- 
12.59); 1 E-04(0-2E-03) 



Data for vaccination are given according to its total coverage percentage and data for treatment, hand-hygiene & clean water distribution (PH & CWD) and sanitation 
given according to average coverage percentage per day. All data are given in the format [mean (95% CI)]. Data for the Bulawayo province is not given as different 
interventions have no effect on case or death reduction in this region. 

Notations in the first column are exactly same as Table 5. * indicate the intervention/intervention combinations which do not have any effect on case or death reduction 
in a province. 

doi:1 0.1 371 /journal.pone.0081 231 .t007 



from improper waste management system [53]. In spite of this; the 
national water supply agency of Zimbabwe (ZINWA) is not in 
charge in supplying water in Bulawayo [54]. 

Mukandavire et. al. [14] estimated R 0 for 2008-2009 cholera 
epidemics in Zimbabwe with constant transmission rate but to our 
knowledge, this is the first time that the basic reproduction 
numbers with periodic transmission rate are estimated for cholera 
epidemics in Zimbabwe or any other country. It is also to be noted 
that our data set contains much longer time scale (from the 
beginning of cholera epidemic in 2008-2009 to June 2011) than 



the earlier analyzed data [14]. It is already pointed out that the 
estimated values of Ro in Bulawayo, Mashonaland West, Manica- 
land, Matabeleland South and Matabeleland North by Mukanda- 
vire et. al. [14] differ with our estimates. The disease dynamics of 
cholera may not be captured properly by assuming constant 
contact rate between human and bacterial populations over time, 
as it also depends on temporal forcing. Thus, the model with 
periodic environment is more appropriate than the previous 
studies. We believe that the prediction, thereby proposed, will be 
helpful for policy makers. 
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Table 8. Optimal cost (in USD) projected between the end of 2008-09 epidemic to January 1, 2012, by province under each 
intervention scenario at an optimal rate. 







Harare 


Mash west 


Mash cen 


Mash east 


Midlands 


PH & CWD 


8.28E+04 (7.39E+04- 
9.56E+04) 


2.82E+05 (2.62E+05- 
3.24E+05) 


3.75E+04 (3.14E+04- 
4.36E+04) 


8.16E+03 (6.28E+03- 
1 .89E+04) 


8.95E+03 (7.43E+03-1 .39 
£+04) 


VA 


3.25E+06 (4.12E+05- 
3.20E+07) 


4.55E+06 (2.28E+06- 
2.34E+07) 


3.20E+06 (5.89E+05- 
1 .85E+07) 


1.91E+06 (1.68E+05- 
1 .87E+07) 


4.07E+06 (1.27E+05- 
3.06E+07) 


TR 


1.26E+05 (1.10E+05- 
1 .48E+05) 


5.13E+05 (4.76E+05- 
5.48E+05) 


2.78E+04 (2.19E+04- 
3.42E+04) 


2.69E+03 (2.17E+03- 
4.40E+03) 


3.07E+03 (2.29E+03- 
3.98E+03) 


SN 


2.96E+05 (2.57E+05- 
3.49E+05) 


1.24E+06 (1.15E+06- 
1 .34E+06) 


8.05E+04 (6.55E+04- 
9.65E+04) 


9.77E+03 (8.34E+03- 
1 .36E+04) 


1.16E+04 (9.49E+03- 
1 .40E+04) 


PH & CWD + TR + VA 


1 .99E+06 (4.48E+04- 
2.67E+07) 


2.22E+06 (1.49E+05- 
2.16E+07) 


1.25E+06 (1.78E+04- 
1 .75E+07) 


3.91E+06 (2.31E+03- 
1 .83E+07) 


3.17E+06 (2.57E+03- 
2.97E+07) 


PH & CWD + TR + SN 


4.97E+04 (4.45E+04- 
5.74E+04) 


1.58E+05 (1.48E+05- 
1 .65E+05) 


2.20E+04 (1.78E+04- 
2.63E+04) 


2.74E+03 (2.31 £+03- 
4.03E+03) 


3.23E+03 (2.41 E+03- 
4.29E+03) 


PH & CWD + VA + SN 


3.99E+06 (7.20E+04- 
3.21 E+07) 


6.50E+06 (2.49E+05- 
2.46E+07) 


1.20E+06 (3.21 £+04- 
1 .76E+07) 


1 .60E+06 (6.32E+03- 
1 .82E+07) 


2.59E+06 (7.99E+03- 
3.12E+07) 


VA + TR + SN 


4.55E+06 (1.02E+05- 
3.29E+07) 


1.09E+06 (4.21 E+05- 
1 .60E+07) 


1 .98E+06 (2.02E+04- 
1 .78E+07) 


9.65E+05 (2.15E+03- 
1.61 E+07) 


4.31 £+06 (2.45E+03- 
3.24E+07) 


PH & CWD + TR+ VA + SN 


1.11E+06 (4.47E+04- 
2.52E+07) 


2.36E+06 (1.48E+05- 
2.11 E+07) 


1.31E+06 (1.78E+04- 
1 .74E+07) 


7.55E+05 (2.31 £+03- 
1.84E+07) 


3.21E+03 (2.41E+03- 
4.34E+03) 




Masvin 


Manica 


Mata south 


Mata north 


Total 


PH & CWD 


6.96E+03 (5.02E+03- 
2.05E+04) 


1.15E+05 (1.06E+05- 
1 .24E+05) 


3.23E+04 (3.03E+04- 
3.43E+04) 


4.26£+04 (3.82E+04- 
4.97E+04) 


6.16E+05 (5.60E+05- 
7.25E+05) 


VA 


2.54E+06 (1.53E+05- 
1 .90E+07) 


6.04E+06 (6.38E+05- 
3.75E+07) 


1 .03E+06 (9.39E+04- 
7.21 £+06) 




2.66E+07 (4.46E+06- 
1 .87E+08) 


TR 


1.53E+03 (1.39E+03- 
1 .72E+03) 


2.03E+05 (1.83E+05- 
2.22E+05) 


3.87E+04 (3.52E+04- 
4.19E+04) 


* 


9.16E+05 (8.31E+05- 
1 .00E+06) 


SN 


6.84E+03 (6.32E+03- 
7.67E+03) 


4.63E+05 (4.18E+05- 
5.10E+05) 


8.29E+04 (7.62E+04- 
9.07E+04) 


* 


2.19E+06 (1.99E+06- 
2.42E+06) 


PH & CWD + TR + VA 


1.61E+06 (1.56E+03- 
1 .72E+07) 


2.66E+06 (6.34E+04- 
2.32E+07) 


1 .36E+06 (2.04E+04- 
7.55E+06) 


2.56E+05 (2.71 £+04- 
2.85E+06) 


1 .84E+07 (3.29E+05- 
1 .64E+08) 


PH & CWD + TR + SN 


1.71 £+03 (1.56E+03- 
1 .87E+03) 


6.77E+04 (6.31E+04- 
7.25E+04) 


2.17E+04 (2.04E+04- 
2.29E+04) 


3.02E+04 (2.71 £+04- 
3.51 £+04) 


3.57E+05 (3.27E+05- 
3.90E+05) 


PH & CWD + VA + SN 


2.23E+06 (5.23E+03- 
1 .75E+07) 


2.87E+06 (1.03E+05- 
2.48E+07) 


1.25E+06 (3.02E+04- 
7.07E+06) 


4.25E+04 (3.81 £+04- 
4.96E+04) 


2.23E+07 (5.44E +05- 
1 .73E+08) 


VA + TR + SN 


1.89E+06 (1.42E+03- 
1 .75E+07) 


3.91E+06 (1.69E+05- 
3.44E+07) 


9.29E+05 (3.45E+04- 
7.21E+06) 


1.86E+05 (4.69E+04- 
1 .72E+06) 


1 .98E+07 (8.00E+05- 
1 .76E+08) 


PH & CWD + TR + VA + SN 


7.20E+05 (1.57E+03- 
1 .76E+07) 


2.98E+06 (6.33E+04- 
2.62E+07) 


2.17E+04 (2.04E+04- 
2.29E+04) 


4.71 £+05 (2.71 £+04- 
1 .04E+07) 


9.74E+06 (3.27E+05- 
1 .36E+08) 



Costs are given in the format [mean(95% CI)]. Cost corresponding to Bulawayo province is not given as different interventions have no effect on case or death reduction 
in this region. 

Notations in the first column are exactly same as Table 5. Here Ek= 10*. * indicate the intervention/intervention combination which do not have any effect on case or 

death reduction in a province. 

doi:1 0.1 371 /joumal.pone.0081 231 .t008 



Optimal cost effective study in Zimbabwe, from the end of 
epidemic in 2008-09 to January 1, 2012, suggests that, as a single 
intervention hand-hygiene & clean water supply is the most cost- 
effective way to control a future cholera outbreak in those regions 
where slow transmission is the dominating factor for cholera 
transmission (see, Tables 8 and 9). In terms of cases and deaths 
reduction during epidemic hand-hygiene & clean water supply is 
by far the most effective individual intervention among the other 
single interventions (see Tables 5 and 6). This result is in good 
agreement with the observations by Andrews and Basu [6], where 
they argued that hand-hygiene & clean water distributions will 
avert more cases and deaths than treatment and vaccination 
during the epidemic in Haiti. Treatment is the most cost-effective 
in those regions where hyper-infectious transmission is the main 



factor for cholera transmission (see, Tables 8 and 9). This result is 
in well agreement with the previous observation of Naficy et. al. 
[55] on the control of cholera in sub-Saharan refugee settings. 
Treatment and hand-hygiene & clean water supply with any other 
intervention combination will also be cost-effective, which could 
avert thousands of cases and hundreds of deaths in Zimbabwe at a 
minimal cost. A synchronized, timely and efficient intervention 
will effectively reduce the severity of the disease and number of 
deaths. Our mathematical model and its prediction will help the 
public health authority of Zimbabwe for making suitable 
intervention strategies. We also believe that, a similar method 
can be applied for endemic and epidemic cholera outbreaks in 
other regions/countries as well, in particular, with seasonal 
patterns in disease transmission. 
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Table 9. Cost per averted case (in USD) projected between the end of 2008-09 epidemic to January 1, 2012, by province under 
each intervention scenario at an optimal rate. 







Harare 


Mash west 


Mash cen 


Mash east 


Midlands 


PH & CWD 


1.30E+02 (1.17E+02- 
1 .42E+02) 


1.01E+02 (9.40E+01- 
1 .23E+02) 


4.02E+01 (3.35E+01- 
4.59E+01) 


5.92E+01 (4.70E+01-1.31E+02) 


5.73E+01 (4.97E+01-9.63 
E+01) 


VA 


5.30E+05 (4.73E+04- 
8.11 £+05) 


1 .97E+05 (2.48E+03- 
6.85E+05) 


1 .79E+04 (1 .99E+03- 
1 .07E+05) 


1.26E+05 (9.88E+04-1.55E+05) 


1.78E+05 (1.17E+05- 
2.19E+05j 


TR 


3.22E+02 (2.89E+02- 
3.52E+02) 


2.62E+02 (2.42E+02- 
2.82E+02) 


3.09E+01 (2.41 E+01 - 
3.68E+01) 


1.95E+01 (1.74E+01-2.76E+01) 


1.95E+01 (1.71E+01- 
2.24E+0/j 


SN 


9.99E+02 (8.96E+02- 
1 .09E+03) 


7.77E+02 (7.21 E+02- 
8.41 E+02) 


9.44E+01 (7.42E+01- 
1.12E+02) 


7.24E+01 (6.66E+01-8.84E+01) 


7.65E+01 (7.01E+01- 
8.54£+0/j 


PH & CWD +TR+VA 


3.12E+03 (6.89E+01- 
4.16E+04) 


7.87E+02 (5.25E+01- 
7.66E+03) 


1.42E+03 (1.91E+01- 
2.06E+04) 


2.80E+04 (1.81E+01-1.35E+05) 


2.16E+04 (1.82E+0/- 
2.)0E+05) 


PH & CWD + TR + SN 


7.70E+01 (6.87E+01- 
8.40E+01) 


5.61E+01 (5.20E+01- 
5.99E+01) 


2.36E+01 (1.91E+01- 
2.74E+01) 


1.98E+01 (1.79E+01-2.53E+01) 


2.04E+01 (1.80E+0/- 
2.6 /E+01) 


PH & CWD + VA + SN 


6.03E+03 (1.15E+02- 
4.96E+04) 


2.25E+03 (8.86E+01- 
8.71E+03) 


1.36E+03 (3.24E+01- 
2.04E+04) 


1.12E+04 (4.55E+01-1.28E+05) 


1.80E+04 (5.19E+0/- 
2.03E+05) 


VA + TR + SN 


8.25E+03 (2.45E+02- 
6.00E+04) 


4.72E+02 (2.00E+02- 
6.39E+03) 


2.23E+03 (2.25E+01- 
2.09E+04) 


6.60E+03 (1.70E+01-1.12E+05) 


2.94E+04 (1.75E+01- 
2.24E+05) 


PH & CWD + TR + 
VA + SN 


1.87E+03 (6.86E+01- 
4.38E+04) 


8.34E+02 (5.20E+01- 
7.53E+03) 


1.50E+03 (1.91E+01- 
2.04E+04) 


5.07E+03 (1.80E+01-1.25E+05) 


2.02E+01 (1.78E+01- 
2.41 E+01) 




Masvin 


Manlca 


Mata south 


Mata north 




PH & CWD 


5.18E+01 (3.31E+01- 
1 .54E+02) 


1.76E+02 (1.67E+02- 
1 .88E+02) 


2.22E+02 (2.06E+02- 
2.34E+02) 


8.18E+02 (7.26E+02-9.06E+02) 




VA 


1.19E+05 (9.24E+04- 
1 .39E+05) 


1 .05E+06 (4.93E+04- 
1 .50E+06) 


1.86E+05 (4.10E+04- 
2.49E+05) 






TR 


1.14E+01 (9.66E+00- 
1.32E+01) 


5.59E+02 (5.19E+02- 
6.02E+02) 


6.29E+02 (5.21 E+02- 
6.95E+02) 






SN 


5.16E+01 (4.42E+01- 
5.93E+01) 


1.72E+03 (1.60E+03- 
1 .86E+03) 


2.00E+03 (1.66E+03- 
2.21 E+03) 






PH & CWD + TR + VA 


1.12E+04 (1.08E+01- 
1 .24E+05) 


3.90E+03 (9.69E+01- 
3.30E+04) 


8.47E+03 (1.40E+02- 
4.92E+04) 


4.21 E+03 (5.45E+02-4.57E+04) 




PH & CWD + TR + SN 


1.27E+01 (1.08E+01- 
1.43E+01) 


1.03E+02 (9.66E+01- 
1.10E+02) 


1.52E+02 (1.40E+02- 
1.61 E+02) 


6.04E+02 (5.37E+02-6.66E+02) 




PH & CWD + VA + SN 


1.56E+04 (3.56E+01- 
1 .25E+05) 


4.23E+03 (1.60E+02- 
3.74E+04) 


7.68E+03 (2.04E+02- 
4.43E+04) 


8.17E+02 (7.26E+02-9.05E+02) 




VA + TR + SN 


1.36E+04 (9.84E+00- 
1 .24E+05) 


7.10E+03 (4.27E+02- 
6.01E+04) 


7.92E+03 (4.67E+02- 
6.02E+04) 


4.56E+04 (1.22E+04-1.12E+05) 




PH & CWD + TR + 
VA + SN 


4.84E+03 (1.08E+01- 
1 .20E+05) 


4.37E+03 (9.66E+01- 
3.85E+04) 


1.52E+02 (1.40E+02- 
1.61 E+02) 


8.01 E+03 (5.46E+02-1.79E+05) 





Data are given in the format [mean (95% CI)]. Data corresponding to Bulawayo province is not given as different interventions have no effect on case or death reduction 
in this region. 

Notations in the first column are exactly same as Table 5. Here Ek= 10*. * indicate the intervention/intervention combination which do not have any effect on case or 

death reduction in a province. 

doi:l 0.1 371 /joumal.pone.0081 231 .t009 
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